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ABSTRACT 

On 7 January 2014 an XI.2 flare and GME with a radial speed ss2500 kms“^ was observed from 
near an active region close to disk center. This led many forecasters to estimate a rapid arrival at 
Earth (si 36 hours) and predict a strong geomagnetic storm. However, only a glancing CME arrival 
was observed at Earth with a transit time of «49 hours and a A"p geomagnetic index of only 3—. We 
study the interplanetary propagation of this CME using the ensemble Wang-Sheeley-Arge (WSA)- 
ENLIL-|-Cone model, that allows a sampling of CME parameter uncertainties. We explore a series 
of simulations to isolate the effects of the background solar wind solution, CME shape, tilt, location, 
size, and speed, and the results are compared with observed in-situ arrivals at Venus, Earth, and 
Mars. Our results show that a tilted ellipsoid CME shape improves the initial real-time prediction to 
better reflect the observed in-situ signatures and the geomagnetic storm strength. CME parameters 
from the Graduated Cylindrical Shell model used as input to WSA-ENLIL-|-Cone, along with a tilted 
ellipsoid cloud shape, improve the arrival-time error by 14.5, 18.7, 23.4 hours for Venus, Earth, and 
Mars respectively. These results highlight that CME orientation and directionality with respect to 
observatories play an important role in understanding the propagation of this CME, and for forecasting 
other glancing CME arrivals. This study also demonstrates the importance of three-dimensional CME 
fitting made possible by multiple viewpoint imaging. 

Keywords: Sun: coronal mass ejections (CMEs), Sun: solar-terrestrial relations, magnetohydrody¬ 
namics (MHD) 


1. INTRODUCTION 

A very fast CME (si 2500 kms“^) on 7 Jan¬ 
uary 2014 was associated with an active region 
near disk center, leading most forecasters world¬ 
wide to predict rapid arrivals at Earth with strong 
geomagnetic impacts on the CME Scoreboard 
(kauai.ccmc.gsfc.nasa.gov/CMEscoreboard). The Com¬ 
munity Coordinated Modeling Center (CCMC) (located 
at NASA Goddard Space Flight Center) serves the CME 
Scoreboard website to the research community who may 
submit CME shock arrival-time predictions in real-time 
using different forecasting methods; this facilitates 
model validation under real-time conditions. Users 
submitted fifteen arrival-time predictions to the CME 
Scoreboard for this event, with an average predicted 
arrival-time of 9 January at 6:35 UT (Ri36 hours transit 
time) and average predicted Kp range of 6-7.6. The 
actual CME arrival-time was at 19:39 UT (Ri49 hours 
transit time) and the geomagnetic Kp index was 3— 
(2.7). 

^ Catholic University of America, Washington, DC, USA 

^ Heliophysics Science Division, NASA Goddard Space Flight 
Center, Greenbelt, MD, USA 

® Department of Astronomy, University of Maryland, College 
Park, MD, USA 

^ Space Science Division, Naval Research Laboratory, Washing¬ 
ton, DC, USA 

® Solar Section, Applied Physics Laboratory Johns Hopkins, 
Laurel, MD, USA 

® George Mason University, Fairfax, VA, USA 

^ Space Research Institute, Austrian Academy of Sciences, Graz, 
Austria 

® IGAM-Kanzelhohe Observatory, Institute of Physics, Univer¬ 
sity of Graz, Graz, Austria 


Mostl et al. (2015) studied the non-radial motion of 
this CME by Ri37°±10° away from the source region 
longitude using EUV and coronagraph observations to¬ 
gether with in-situ signatures at Earth and Mars. They 
demonstrate that the CME was “channeled” by strong 
nearby active region magnetic fields and open coronal 
fields into a non-radial propagation direction within 2.1 
solar radii (R©), in contrast to deflection in interplane¬ 
tary space. Torok & Kliem (2005) and Liu (2008) show 
that instabilities in a modeled magnetic flux rope can 
lead to a CME or a failed eruption depending on how 
the overlying coronal field decreases with height. In a 
recent study by Thalmann et al. (2015) it was shown, 
using reconstructions of the three-dimensional (3D) coro¬ 
nal magnetic field from a nonlinear force free model, that 
the overlying magnetic field plays an important role in 
producing CMEs in the lower corona. Together with 
the result from Mostl et al. (2015), this indicates that 
the coronal magnetic field is an important parameter for 
early CME characteristics. 

In this work we consider the relative importance of 
CME geometry with respect to observatories in CME 
propagation modeling. Whether a CME will arrive at 
Earth, its arrival-time, and geo effectiveness is largely de¬ 
termined by the angle between the CME central axis with 
respect to Earth and the CME size. This is complicated 
by the large uncertainties in CME size, and the related 
properties of speed and direction. Many CME propaga¬ 
tion models do not take into account the full 3D CME 
geometry, such as tilt, asymmetric size, or propagation 
direction out of the ecliptic plane (2D models). Shen 
et al. (2014) found that for the front-sided full halo CMEs 
which arrived at Earth in their sample, 74% of them 
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satisfied the criterion that the CME half angular width 
was greater than the angle between the CME propaga¬ 
tion direction and Sun-Earth line. Most geoeffective halo 
CMEs originate from near the central meridian, but only 
about half of Earth-directed halo CMEs are geoeffective 
(Schwenn et al. 2005; Wang et al. 2002; Kim et al. 2008). 
Moon et al. (2005) show that for a sub-selected sam¬ 
ple of twelve front-sided halo CMEs the CME direction 
is an important parameter for determining the geoeffec¬ 
tiveness, as measured by the Dst index. In their study, a 
“direction parameter”, determined from how asymmet¬ 
ric the halo CME appeared in coronagraph observations 
was well-correlated with Dst, while the source region lo¬ 
cation and CME speed, which is highly affected by pro¬ 
jection effects, were not as well-correlated. Kim et al. 
(2008) confirmed this result for a larger sample of 486 
front-sided halo CMEs. 

Jian et al. (2006) categorized CME total pressure pro¬ 
files can into three regimes depending on the impact an¬ 
gle with respect to the CME core; the impact param¬ 
eter affects the observed magnetic field structures and 
any geomagnetic response. Savani et al. (2015) present 
a method to predict the CME’s magnetic field at Earth 
based on the Bothmer & Schwenn (1994, 1998) helicity 
rule and the impact parameter relative to the CME flux- 
rope axis. They highlight the importance of triangulat¬ 
ing CME geometry when predicting the resulting in-situ 
properties at Earth. 

Uncertainties in CME properties can be reduced, but 
not eliminated, by CME reconstruction techniques which 
take advantage of multiple viewpoint imaging. Mierla 
et al. (2010) evaluate different techniques to reconstruct 
the 3D configuration of CMEs from coronagraph data 
from multiple viewpoints and estimate that the CME 
propagation direction can be determined to within 10°. 
Colaninno et al. (2013) show that CME arrival-time pre¬ 
dictions based on 3D height-time measurements of So¬ 
lar TErrestrial RElations Observatory [STEREO: Kaiser 
et al. 2008) Sun Earth Connection Coronal and Helio¬ 
spheric Investigation (SECCHI: Howard et al. 2008) he¬ 
liospheric images show a half-day improvement compared 
with methods using one coronagraph viewpoint. Mdstl 
et al. (2014) also find that predicting CME speeds and 
arrival-times using techniques applied to J-maps con¬ 
structed from heliospheric images yielded more accu¬ 
rate results than using plane-of-sky coronagraph mea¬ 
surements. Shi et al. (2015) emphasize the importance 
of using stereoscopy to determine CME parameters for 
the accuracy of arrival-time for 21 Earth-directed CMEs 
with the Drag Based Model (DBM: Vrsnak et al. (2004, 
2013)) model. 

In order to further understand the interplanetary prop¬ 
agation of the 7 January 2014 CME we use the en¬ 
semble WSA-ENLIL-|-Cone model implemented at the 
CCMC. The Wang-Sheeley-Arge (WSA) coronal model 
(Arge & Pizzo 2000; Arge et al. 2004) is used as input 
to the ENLIL (Odstrcil 2003) 3D magnetohydrodynamic 
(MHD) code that computes the background solar wind 
into which a hydrodynamic “CME” cloud is launched. 
Ensemble modeling is a collection of forecasts for the 
same event that represent possible scenarios given the 
uncertainties associated with the forecasting (such as ini¬ 
tial conditions, techniques, and models). In addition to 
giving an average forecast, probabilities of events can be 


determined from the possible scenarios. In this study we 
evaluate the sensitivity of the WSA-ENLIL-fCone mod¬ 
eled CME arrival-times to the uncertainties associated 
with the initial CME parameters of location, speed, size, 
and shape. 

This paper is organized as follows: In Section 2 we 
discuss the observations and methods for analyzing CME 
kinematics to be used as input to the model. The WSA- 
ENLIL-|-Cone ensemble model is described in Section 3. 
A series of ensemble CME propagation modeling results 
are presented in Sections 3.1-3.2, and further compared 
in Section 3.3. Finally, we present a summary of the 
results in Section 4. 

2. OBSERVATIONS 

The 7 January 2014 CME was associated with an XI.2 
class solar flare from between active regions numbered 
11943 and 11944 located at S12°W08° and peaked at 
18:32 UT. Figure la shows the Solar Dynamics Obser¬ 
vatory [SDO: Pesnell et al. 2012) Atmospheric Imag¬ 
ing Assembly (AIA: Lemen et al. 2012) 193 A difference 
image at 19:00 UT from which the 18:00 UT image is 
subtracted (prior to flare). The SDO/Helioseismic and 
Magnetic Imager (HMI: Scherrer et al. 2012) colorized 
magnetogram image is superimposed on this difference 
image. In the colorized magnetogram overlay, yellow-red 
colors indicate negative polarity, and green-blue colors 
indicate positive polarity. Bright post-eruption loops are 
seen near the active region with an extensive dimming 
region visible to the south and southwest, with weak neg¬ 
ative polarity. Figure lb shows the SDO/AIA 211A im¬ 
age near 18:00 UT with open (pink: negative polarity, 
green: positive polarity) and closed (white) magnetic 
field lines from a Potential Field Source Surface (PFSS) 
model overlaid (Schrijver et al. 2006). Two large coronal 
holes are visible to the northeast of the active region, and 
open coronal field is present to the west of the region. 

The top row of Figure 2 shows the CME first ob¬ 
served at 18:24 UT on 7 January 2014 by the SOlar 
and Heliospheric Observatory [SOHO: Domingo et al. 
1995) Large Angle and Spectrometric Coronagraph Ex¬ 
periment (LASCO: Brueckner et al. 1995) C3 and 
STEREO/SDCCDl COR2-A and -B coronagraphs near 
the time of 19:40 UT. STEREO-A was located at 151° 
longitude (west of the Sun-Earth line), and STEREO B 
at —153° longitude (east of the Sun-Earth line), in He¬ 
liocentric Earth Equatorial (HEEQ) coordinates®. Earth 
was located at 0° longitude, —3.8° latitude; Venus at 
—2.6° longitude, —1.8° latitude; and Mars at 52° longi¬ 
tude, —5.4° latitude. 

Three other energetic CMEs occurred in the days lead¬ 
ing up to 7 January. An Earth-directed CME originat¬ 
ing from near active region number 11943 was first ob¬ 
served on 4 January 2014 at 19:39 UT by STEREO- 
A/SECCHI COR2 and was associated with an M4.0 
flare from S09°E34° peaking at 19:46 UT. The CME 
was propagating Ri 4° west of the Sun-Earth line, —35° 

® In this coordinate system the Z axis is aligned with the solar 
north rotation pole and the X axis pointing toward the intersection 
between the solar equator and the solar central meridian as seen 
from Earth (Hapgood 1992; Thompson 2006). HEEQ coordinates 
are related to Stonyhurst heliographic coordinates, with directions 
south of the origin represented by negative HEEQ latitudes, and 
directions east by negative HEEQ longitudes. 
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Figure 1. (a) SDO/HMI colorized magnetogram image superimposed on a SDO/AIA 193 A difference image at 2014-01-07 19:00 UT 

from which the 18:00 UT image is subtracted (prior to flare). In the colorized magnetogram overlay, yellow-red colors indicate negative 
polarity, and green-blue colors indicate positive polarity. Bright post-eruption loops were seen near the active region with an extensive 
dimming region visible to the south and southwest of the active region with weak negative polarity, (b) SDO/AIA 21lA image near 18:00 
UT with open (pink: negative polarity, green: positive polarity) and closed (white) magnetic field lines from a PFSS model overlaid. Two 
large coronal holes are visible to the northeast of the active region, and open coronal field is present to the west of the region. 



Figure 2. Top row: Coronagraph observations of the 7 January 2014 CME as viewed from (a) STEREO-B SECCHI/COR2, (b) SOHO 
LASCO/C3, and (c) STEREO-A SECCHI/COR2, near the time of 19:40 UT. The fields-of-view of LASCO C3 and SECCHI/COR2 
are 2.8-32 R© and 2.5-15 Rq respectively. Bottom row: Observations from the same instruments with the Time Convolution Mapping 
Method applied from 18:24-20:39 UT. Each color represents a different STEREO-A SECCHI/COR.2 image time from 18:24-20:39 UT, 
with STEREO-B SECCHI/COR2 and SOHO LASCO/C3 images synchronized within 6 minutes of the STEREO-A time. The solid arrows 
point out sheath/shock-like structures identified with TCMM, and dotted line arrows point out other parts of the CME that show different 
kinematic profiles from the “driver” or “main-body” such as along the CME legs. 
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south of the solar equatorial plane, with a speed of 
R::!750±200 kms“^. Later, on 7 January at 14:24 UT 
this CME arrived at Earth but did not produce a ge¬ 
omagnetic storm. Shortly after, on 4 January at 23:08 
UT {STEREO-A/SECCRl COR2) a CME erupted from 
active region 11936 just beyond the west limb with a 
speed of w 650 ±100 kms“^ associated with a partially 
occulted Ml.9 flare. A small increase in the background 
solar energetic particle flux was observed near Earth 
by the Geostationary Operational Environment Satel¬ 
lites) (GOES)-13/Electron, Proton, Alpha, and Detector 
(EPEAD) starting on 5 January. Another CME origi¬ 
nating from just beyond the west limb was observed on 

6 January 2014 at 08:00 UT by 50770/LASC0 C2 from 
the same active region (11936), with an average speed 
of ^llOO kms“^. The CME reached a peak speed of 
«1960 kms“^ and was associated with a partially oc¬ 
culted C2.2 flare, a solar energetic particle event, and 
a ground level enhancement (Thakur et al. 2014). The 
solar energetic particle intensity remained enhanced be¬ 
fore increasing again due to the XI.2 solar flare and 

7 January CME, to above Ipfu/MeV (1 pfu = l parti¬ 
cle cm“^ sr“^ s“^) as observed near Earth by the GOES- 
IJ'/EPEAD 15-40MeV proton channel. 

2.1. Time Convolution Mapping Method 

In this study, the CME ejecta (driver) is the input to 
the WSA-ENLIL±Cone model, therefore it is important 
to separate the CME ejecta from CME-associated bright- 
enings such as streamer deflections and compressive wave 
fronts when interpreting the coronagraph images. We ap¬ 
plied a CME identification algorithm. Time Convolution 
Mapping Method (TCMM; Thompson & Young (submit¬ 
ted, 2015); Thompson et al. (2015)), to the coronagraph 
data (top row of Eigure 2) to identify the time history 
of the CME evolution. The advantage to TCMM is that 
the spatiotemporal evolution of a CME can be captured, 
allowing users to separate features with different prop¬ 
agation characteristics. For example, separating “true” 
or “main” CME mass from CME-associated brightenings 
is a well-known obstacle (Gopalswamy et al. 2009; Kwon 
et al. 2015). A TCMM CME map is made by first record¬ 
ing the maximum value each individual pixel in the im¬ 
age reaches during the traversal of the CME (in this case 
from 18:24 UT to 20:39 UT). Then the maximum value 
is convolved with a color to indicate the time that the 
pixel reached that value. The lower row of Figure 2 shows 
the TCMM maps corresponding to the viewpoints in the 
top row for multiple time intervals during the 7 January 
2014 CME propagation. The TCMM user is then able to 
identify continuous “rainbow profiles,” indicating related 
kinematic behavior, and also identify breaks in the rain¬ 
bow profiles that indicate a discontinuity in kinematic 
history. Each color represents a different STEREO-A 
SECCH1/COR2 image time from 18:24-20:39 UT, with 
STEREO-B SECCm/com and 50770 LASCO/C3 im¬ 
ages synchronized within 6 minutes of the STEREO time. 
Continuous loop-like TCMM structure (most visible in 
the lower right STEREO-A frame of Figure 2) is indica¬ 
tive of the CME driver as described by Gopalswamy et al. 
(2009). Kinematic rainbow profiles that appear as con¬ 
tinuous loops are in contrast to the profiles with discon¬ 
tinuities that are more consistent with radial structures 


seen outside of the CME ejecta. These discontinuities are 
more indicative of the interaction of the erupting struc¬ 
ture with the surrounding corona. This method also dis¬ 
tinguishes the faint outer structure ahead of the bright 
CME loop, which some identify as the CME-driven shock 
front or compressive wave (Vourlidas et al. 2003, 2013; 
Savani et al. 2013). Kwon et al. (2015) have found evi¬ 
dence that many CMEs appear as full halos because of 
this surrounding structure, and not due to the viewing 
geometry of the main CME driver. The solid arrows in 
Figure 2 point out sheath/shock-like structures identi¬ 
fied with TCMM, and dotted line arrows point out other 
parts of the CME that show different kinematic profiles 
from the “main body”, such as along the CME legs. The 
STEREO-A TCMM map also shows how the northern 
portion of the CME is deformed such that its motion is 
impeded compared with the southern CME front. 

2.2. CME Kinematics 

Many different techniques exist for deriving the CME 
parameters of location, speed, and size that are used as 
input to the WSA-ENLIL±Cone model. In this section 
we consider a variety of techniques to create an input 
parameter distribution that adequately samples the error 
bar range for this CME. Each distribution is then used 
as input to the ensemble simulations A and B which will 
be described in Section 3. 

CCMC derived a first round of true 3D parameters 
(de-projected) in real-time for the 7 January CME, 
using plane-of-sky measurements of beacon data from 
50i70/LASC0-C3, STEREO/SECCBl COR2 com¬ 
bined with the CME direction that was estimated from 
the extensive dimming to the southwest of the active re¬ 
gion. The de-projected speed was calculated using the ra¬ 
tio of projected speed to the true speed relation from the 
theoretical geometric model by Hundhausen et al. (1994) 
using the CME direction and projected CME half-width 
(57°-72°). With this technique the CME propagation 
direction was measured to range from —23° to —31° lat¬ 
itude, 22° to 50° west of the Sun-Earth line in longitude, 
and speed from Ri 2000 to 3000 kms“^. Median CME 
parameters are: speed of 2400 kms“^, direction of 38° 
longitude, —28° latitude, and a half-width of 64°. These 
kinematic parameters are extrapolated for a height of 
21.5 R©, the inner boundary of the ENLIL model. Dis¬ 
crete samples from the range of measured CME param¬ 
eters are shown in Figure 3 in the (a) equatorial plane 
(latitude = 0°) and (b) meridional plane (longitude = 0°) 
and is used as input to ensemble A. The plots show the 
CME velocity vectors in spherical HEEQ coordinates, 
the grids show degrees longitude (a) and latitude (b), 
and the radial coordinate shows the speed in kms“^. 
The Sun-Earth line is along 0° longitude and —3.8° lat¬ 
itude. Arrow directions indicate the CME central lon¬ 
gitude and latitude, with the CME half-width indicated 
by the color of the vector, and arrow lengths correspond 
to CME speed. 

Mdstl et al. (2015) discuss the 3D kinematics of the 
7 January CME using science level data with a vari¬ 
ety of analysis methods. The techniques and param¬ 
eters are summarized in Table 1. Measurements be¬ 
tween 2.1 and 18.5 R© using the Graduated Cylindri¬ 
cal Shell (GCS) model from Thernisien et al. (2006) 
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Table 1 

Summary of de-projected CME parameters from various techniques described in Section 2.2. 


Technique 

V 

(kms“^) 

Lat 

(°) 

Long 

(“) 

a/2 

(“) 

fi/2 

n 

Tilt 

(“) 

Real-time (Ensemble A): 

Hundhausen et al. (1994) geometric model 

2400+ 

-281? 

38 lie 

64 ± 8 

- 

- 

Science level data: 

GCS 

2565 ± 250 

-25 ± 5 

32 ± 10 

35 ± 4 

22 ± 4 

38 ± 5 

Fixed-0 and Harmonic Mean* (triangulation) 

2124 ± 283 

- 

45 ± 10 

- 

- 

- 

Spherical Shell 

2000 ± 200 

-28 ± 2 

32 ± 2 

- 

- 

- 

Combination of methods (Ensemble B): 


oicy + 

340 

-25 ± 5 

36 ± 14 

44 ± 12 

28 ± 8 

38 


* parameters are for the ecliptic component of the CME. 




speed (km/s) .go* 


Median parameters: 
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Figure 3. Ensemble A: Distribution of the 7 January 2014 CME kinematic properties derived in real-time (using plane-of-sky estimates 
combined with the dimming location) shown in the (a) equatorial plane (latitude =0 °) and (b) meridional plane (longitude = 0°). The 
plots show the CME speed vectors in spherical HEEQ coordinates with the grids showing the degrees longitude (a) and latitude (b), and 
the radial coordinate showing the speed in kms“^. The Sun-Earth line is along 0° longitude and —3.8° latitude and the Sun-Mars line 
along 52° longitude and —5.4° latitude. Arrow directions indicate the CME central longitude and latitude, with the CME half-width 
indicated by the color of the vector, and arrow lengths correspond to CME speed. The CME propagation direction ranges from —23° to 
—31° latitude, 22° to 50° west of the Sun-Earth line in longitude, and speed ranges from Pi 2000 to 3000 kms“^. Median CME parameters 
are: speed of 2400 kms“^, direction of 38° longitude, —28° latitude, and a half-width of 64°. 
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Figure 4. Ensemble B: Distribution of the 7 January 2014 CME kinematic properties derived from the CCS model and other analysis 
techniques, in the same format as Figure 3. The CME ellipse aspect ratio is 1.59 with the major axis half-width indicated by the color, 
and CME tilt is 38° counter clock-wise relative to the solar equator. The CME propagation direction ranges from —20° to —30° latitude, 
22° to 49° west of the Sun-Earth line in longitude, and speed ranges from 1817 to 2685 kms“^. Median CME parameters are: speed of 
2157 kms~^, direction of 36° longitude, —25° latitude, major half-width of 44°, and minor half-width of 28°. 
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yielded a speed of 2565 ±250 kms“^, 32° ± 10° longi¬ 
tude, —25°±5° latitude, a/2 = 35° ±4° major axis half- 
width, /3/2 =22° ± 4° minor axis half-width, with a tilt of 
38° ± 5° counter clock-wise relative to the solar equator. 
(The error bars are based on multiple GCS fits and are 
thus particular to this CME configuration). The GCS 
model is a forward modeling technique which assumes 
a magnetic flux rope topology (Thernisien et al. 2006, 
2009). The ratio of the CME face-on width (ellipse ma¬ 
jor axis) to the edge-on width (ellipse minor axis) gives an 
ellipse aspect ratio of 1.59 ± 0.2. Mostl et al. (2015) used 
the in-situ arrival-times at Earth and Mars as constraints 
in the Ellipse Evolution model (ElEvo) model and found 
the ellipse aspect ratio to be in the range of 1.4 ±0.4 
(with no tilt). The average results of the fixed-(/ and Har¬ 
monic Mean triangulation techniques (Liu et al. 2010) be¬ 
tween 20-30 R© give a speed of 2124±283 kms“^, and 
longitude of 45° ± 10° for the ecliptic component of the 
CME (Mostl et al. 2015). We also fit a spherical shell 
model to this CME, and found a speed of 2000 ± 200 
kms“^, 32° ± 2° longitude and —28° ±2° latitude. 

Figure 4 (same format as Figure 3) shows the distribu¬ 
tion of the 7 January CME kinematic properties at 21.5 
Rq from the analysis techniques described above and is 
used as input to ensemble B. This distribution samples 
the error bar range of the CME parameter space derived 
predominantly from the GCS technique, but also from 
HI analysis techniques (see Table 1). The CME ellipse 
aspect ratio is 1.59 with the major axis half-width indi¬ 
cated by the color, and CME tilt is 38° counter clock¬ 
wise relative to the solar equator. The CME propagation 
direction ranges from —20° to —30° latitude, 22° to 49° 
west of the Sun-Earth line in longitude, and speed ranges 
from 1817 to 2685 kms“^. The major and minor half¬ 
widths range from 22° to 54° and 13° to 32° respectively. 
Median CME parameters are: speed of 2157 kms“^, di¬ 
rection of 36° longitude, —25° latitude, major half-width 
of 44°, and minor half-width of 28°. 

3. INTERPLANETARY PROPAGATION 

To study the interplanetary propagation and arrival¬ 
time of this CME we use the ensemble WSA- 
ENLIL±Cone model from the CCMC. The global 3D 
MHD ENLIL model provides a time-dependent descrip¬ 
tion of the background solar wind plasma and magnetic 
field using the WSA coronal model (Arge & Pizzo 2000; 
Arge et al. 2004) as input at the inner boundary of 
21.5 Rq (Odstrcil et al. 1996; Odstrcil & Pizzo 1999a,b; 
Odstrcil 2003; Odstrcil et al. 2004). A homogeneous, 
over-pressured hydrodynamic plasma cloud is launched 
through the inner boundary of the heliospheric compu¬ 
tational domain and into the background solar wind. The 
modeled CME cloud is approximated by a sphere, how¬ 
ever ENLIL also supports other CME shapes such as an 
ellipsoid that can have an arbitrary tilt with respect to 
the solar equator. These shapes can also be elongated 
in the radial direction, and different leading and trail¬ 
ing edge velocities are also possible. In this study, a 
spherical CME cloud (default) and tilted ellipsoid cloud 
are used. Cloud parameters of location, speed and size 
can be determined using any technique for deriving CME 
kinematic properties from coronagraph data. A common 
technique is to assume that the geometrical CME prop¬ 
erties are approximated by the Cone model (Zhao et al. 


2002; Xie et al. 2004) which assumes isotropic expan¬ 
sion, radial propagation, and constant CME cone angular 
width. 

Ensemble modeling of a CME is a collection of CME 
arrival-time predictions that represent possible scenar¬ 
ios given the uncertainties associated with initial condi¬ 
tions from observations and modeling (Lee et al. 2013; 
Emmons et al. 2013; Mays et al. 2015). In addition 
to giving an average CME arrival-time prediction, this 
probabilistic forecast allows one to derive an arrival-time 
uncertainty from the spread and distribution of predic¬ 
tions. The probability of CME arrival can be derived 
from the number of CME arrival predictions as a per¬ 
centage of ensemble size. The current ensemble model¬ 
ing implementation at CCMC evaluates the sensitivity 
of WSA-ENLIL±Cone model simulations of CME prop¬ 
agation to initial CME parameters. For n initial CME 
parameters, n simulations are performed that provide n 
profiles of MHD quantities (density, velocity, tempera¬ 
ture, and magnetic field) and a distribution of n model 
arrival-times at locations of interest within the compu¬ 
tational domain. For each simulation and location of 
interest, a CME arrival “hit” is defined by an increase in 
the dynamic pressure. The ratio Uhits/^ gives a forecast 
probability and conveys the forecast uncertainty about 
whether or not a CME will arrive at that location. For 
this study n = 48 ensemble members are used such that 
48 individual simulations compose each ensemble. De¬ 
tails of the WSA-ENLIL±Cone ensemble modeling sys¬ 
tem at the CCMC are described in Mays et al. (2015). 

For Earth-directed CMEs, n estimates of the ge¬ 
omagnetic ATp index are computed using the WSA- 
ENLIL±Cone model plasma parameters at Earth with 
a modihed Newell et al. (2007) coupling function (Savani 
et al. submitted, 2015; Mays et al. 2015). The function 
represents the rate of magnetic flux (d<i>Mp/dt) opening 

at the magnetopause d<I>Mp/dt = Wbuik^^^7?T'^^sin®/^(^) 
where Ubuik is the bulk solar wind speed, the interplan¬ 
etary magnetic field (IMF) clock angle Oq is given by 
{By/Bz)I and the perpendicular component of the 
magnetic field is given by Rt = {By ± (in GSM 

coordinates). An exponential fit to the correlation of this 
coupling function with the iLp index yields the relation 
K'p = 9.5 — 

We use the 7 January CME parameter distributions 
described in Section 2.2 as input to WSA-ENLIL±Cone 
ensemble simulations and compare these with observa¬ 
tions in the following Sections 3.1-3.2. Ensemble A uses 
the default spherical CME cloud shape and the CME 
parameter distribution derived in real-time using the 
Hundhausen et al. (1994) geometric model (see Figure 

3 and Table 1). Ensemble B, in contrast uses an ellip¬ 
soid CME shape and an input parameter distribution 
predominantly based on the GCS technique (see Figure 

4 and Table 1). Both ensembles A and B use the same 
input WSA background solar wind map. 

3.1. Ensemble A: Real-time Measurements 

CCMC completed the first ensemble simulation-A-for 
this CME in real-time (prior to CME arrival) using the 
input distribution shown in Figure 3 as initial condi¬ 
tions for 48 WSA-ENLIL±Cone simulations. A spherical 
CME shape was used, representing a symmetric CME an- 
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Figure 5. Snapshot of the ENLIL model input inner boundary map at 21.5 Rq computed from the WSA model for the (a) radial velocity 
and (b) number density used as input for one member of ensemble A (left) and ensemble B (right). The map shows the time when the 
inserted CME cross-section is largest, for the (left) spherical CME shape (equal to the input CME half-width) for ensemble A and (right) 
ellipsoid CME shape (equal to the input CME minor and major axis half-widths) from CCS measurements used in ensemble B. Animations 
available in online version. 


gular width. Ensemble A uses the default settings avail¬ 
able from CCMC’s Runs on Request. The input WSA 
model map for this event was computed from a single 
National Solar Observatory Global Oscillation Network 
Group (GONG: Harvey et al. 1996) daily-updated syn¬ 
optic magnetogram for Carrington rotation number 2135 
and Carrington longitude 40° on 7 January 2014 at 18:04 
UT (Earth Carrington longitude was 100°). Figure 5 
(left) shows a snapshot of the ENLIL model input in¬ 
ner boundary map at 21.5 Rq computed from the WSA 
model for the (a) radial velocity and (b) number density 
in HEEQ coordinates, for ensemble A’s median CME 
input parameters. The map shows the time when the 
inserted CME spheroid cross-section, a circle, is largest 
(equal to the input CME half-width) in the background 
solar wind. Red or blue outlines represent outward or in¬ 
ward IMF polarity respectively, and the white line shows 
the heliospheric current sheet (HCS). Figure 6 shows 
a radial velocity contour plot using this input for the 
(a) constant Earth HEEQ latitude plane, (b) meridional 
plane of Earth, and (c) 1 AU sphere in cylindrical projec¬ 
tion on 9 January at 06:00 UT. Panel (d) shows the mea¬ 
sured OMNI hourly average (red) and simulated (blue) 
radial velocity profiles at Earth, with the simulated CME 
duration shown in yellow. The figure shows the north¬ 
eastern portion of the CME impacting Earth. An ani¬ 
mation of the input inner boundary map, and resulting 
simulation are available as electronic supplementary ma¬ 
terial. 


Figures 7a, 8a, and 9a show the modeled magnetic 
field, radial velocity, density, and temperature profiles 
at Venus, Earth, and Mars plotted as color traces for all 
48 ensemble members. The model traces are color-coded 
by CME input speed; slow to faster speeds range from 
light green to dark blue. For Earth the in-situ Li ob¬ 
servations from Advanced Composition Explorer {ACE: 
Stone et al. 1998) are plotted in black (red for B^) and 
Wind is plotted in grey. At Venus, the in-situ observa¬ 
tions from Venus Express ( VEX: Svedhem et al. 2009) 
are shown for periods when the spacecraft was in the so¬ 
lar wind. For UFA, the total magnetic field is shown in 
black and in red. The range of model arrival-times 
for hits (defined by an increase in the model dynamic 
pressure) is indicated by the grey bar on the time axis of 
each figure. The results show that the probability that 
the CME will arrive at Venus, Earth, and Mars is 100% 
(48 hits out of 48 members). 

Figures 7a-9a show that all ensemble A members ar¬ 
rive prior to the observed time at the three locations, 
indicated by the vertical dashed red line. First, VEX 
observed the beginning of the CME-associated distur¬ 
bance between 4:30-7:01 UT (5:45 UT ± 1.25 hours) on 
9 January, as shown in Figure 7a. The VEX data be¬ 
tween about 9 January 4:32 UT (data gap follows) and 
6:35 UT is during a Venus close encounter. After this 
period, a shock is observed with a prompt rise in mag¬ 
netic field to about 19 iiT at 7:01 UT. There is no plasma 
data available at, and about 4 hours around the shock. 
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Figure 6. Ensemble A global view of the 7 January 2014 CME on 9 January at 06:00 UT: WSA-ENLILJ-Cone radial velocity contour 
plot for the (a) constant Earth latitude plane, (b) meridional plane of Earth, and (c) 1 AU sphere in cylindrical projection, for the ensemble 
member with median CME input parameters (speed of 2400 kms“^, direction of 38° longitude, —28° latitude, and a half-width of 64°). 
Panel (d) shows the measured (red) and simulated (blue) radial velocity profiles at Earth, with the simulated CME duration shown in 
yellow. Animation available in online version. 


After a sheath region, there are some magnetic field rota¬ 
tions observed from approximately 9 January 19:00 UT 
to 10 January 10:00 UT. However, the rotations are on 
a small angular scale and there is no well-defined flux 
rope or magnetic cloud observed. At Venus, the ENLIL 
ensemble mean arrival-time was on 8 January at 14:52 
UT with individual arrival-times ranging from -|-7/ — 5 
hours (indicated by the grey bar in Figure 7a) around 
the mean. Using the mean arrival-time prediction, the 
prediction error is Atei-I- — tpi-edlcted ^observed — 14.9 

hours (early prediction). Figure 7a shows that the mod¬ 
eled radial speed is a factor of ~ 1.5 to 2.3 times higher 
than observed by VEX. 

Next, the arrival of the CME was observed by Wind 
and ACE at LI ahead Earth on 9 January at around 
19:39 UT as a weak interplanetary shock (see Figure 
8a). The sheath region after the shock lasts until about 
10 January 07:00 UT, during which the alpha particle 
abundance is increased. The magnetic obstacle is ob¬ 
served from approximately 10 January 07:00 UT to 11 
January 03:00 UT, indicated by magnetic field rotations 
and bidirectional suprathermal electron flux (not shown). 
The magnetic field rotation directions indicated by the 
changes of cone and clock angles are roughly consistent 
with those observed at VEX (not shown). The rotations 
are on a larger scale at Earth than at VEX, however they 
are not smooth enough to be a magnetic cloud (lower- 
than-expected /3 and proton temperature are also not 
observed). The mean ENLIL arrival-time at Earth was 
on 9 January at 00:10 UT with individual arrival-times 
ranging from 4-9/ — 7 hours (indicated by the grey bar 
in Figure 8a) around the mean. Using the mean arrival¬ 
time prediction, the prediction error is —19.4 hours (early 
prediction). Figure 8a shows that the modeled peak solar 
wind speed at the CME arrival-time ranged from 600 to 
1200 kms“^ while the observed value was Ri 500 kms“^. 


The density was also overpredicted by about a factor of 
3 to 5. The model also does not capture the high speed 
stream on 12 January reaching ri 850 kms“^ from the 
coronal holes to the northeast of the active region (see 
Figure lb). The WSA-ENLIL background solar wind 
model driven by GONG synoptic magnetogram shows a 
slower ft! 450 kms“^ southern extension of the coronal 
hole crossing Earth on 14 January. Figure 8a also shows 
the observed arrival of an earlier CME from 4 January 
(described in Section 2) at around 14:24 UT on 7 Jan¬ 
uary. 

The maximum Kp index was calculated using the 
ENLIL modeled quantities of the arrival at Earth (col¬ 
ored traces shown in 8a). Figure 10a shows the pre¬ 
dicted Kp probability distribution for three IMF clock 
angle scenarios 9c = 90° (green), 135° (purple), 180° 
(orange). The figure also shows the overall predicted 
Kp probability distribution calculated for all three an¬ 
gles combined 90°-180°, assuming each equal weights, in 
black. The predicted Kp probability distribution has an 
overall standard deviation cr of 1.8 and there is a ft 78% 
chance for the maximum Kp to be between 6 and 8. Us¬ 
ing the average Kp prediction of 7 or the most likely Kp 
prediction of 8, the Kp prediction error for this event is 
XKp err — Kp predicted Kp observed — 4.33 Or 5.33 (over- 
prediction). For this ensemble simulation, the overpre¬ 
diction of the CME speed and density at Earth leads to 
an overprediction of Kp. Mays et al. (2015) also dis¬ 
cuss a general bias for the overprediction of Kp with this 
method for CME input speeds above ft 1000 kms“^. 

Next, a Forbush decrease was observed at Mars by the 
Mars Science Laboratory (MSL) Radiation experiment 
(RAD) at around 22:30 UT ± 1 hour on 10 January 
and the CME arrival is thought to be within 2 hours of 
the Forbush decrease time (as detailed by Mostl et al. 
(2015)). At Mars, the ENLIL ensemble mean arrival- 
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(a) WSA-ENLIL+Cone Ensemble A Profiles and VEX Observations at Venus (b) WSA-ENLIL+Cone Ensemble B Profiles and VEX Observations at Venus 




Figure 7. Modeled magnetic field, radial velocity, density, and temperature profiles at Venus for (a) ensemble A and (b) ensemble B 
plotted as color traces (slow to fast CME input speeds range from light green to dark blue) for all 48 ensemble members. In-situ observations 
from VEX are plotted black for periods when the spacecraft was in the solar wind. The total magnetic field is shown in black and Bz is 
in red. The vertical dashed red line indicates that the beginning of the CME-associated disturbance is between 4:30-7:01 UT (5:45 UT di 
1.25 hours) on 9 January. The range of modeled arrival-times for hits is indicated by the grey bar on the time axis. 


time was on 9 January at 18:22 UT with individual 
arrival-times ranging from -1-8/ — 8 hours around the 
mean. Using the ensemble mean arrival-time, the pre¬ 
diction error is —28.1 hours (early prediction). Figures 6 
and 9a also shows a high speed stream modeled to arrive 
mid-day on 8 January, just ahead of the modeled CME 
arrival. Ensemble A modeled arrival-times are summa¬ 
rized in Table 2. 

3.2. Ensemble B: GCS Measurements 

Since the GCS model provides the most reliable CME 
3D reconstruction, it is anticipated that modeling CME 
propagation using GCS-derived parmeters, including tilt, 
should be more accurate when comparing to observa¬ 
tions. To explore this, we completed another ensemble- 
B-simulation for this CME with parameters derived from 
GCS measurements and other analysis techniques de¬ 
scribed in Section 2.2 and shown in Figure 4, using an 
ellipsoid CME shape with a tilt. This (and subsequent 
ensembles) use the same input WSA model synoptic map 
as Ensemble A. Figure 5 (right) shows a snapshot of the 
ENLIL model input inner boundary map for a time when 
the inserted tilted CME ellipsoid cross-section (ellipse) is 
largest (equal to the input CME minor and major axis 
half-widths), for one ensemble B member predicting a 
“hit”. The ensemble B member shown has a speed of 
2061 kms“^, direction of 30° longitude, —24° latitude, 
49° major axis half-width, and 31° minor axis half-width. 
For this model input. Figure 11 shows a radial velocity 
contour plot in the same format as Figure 6. This sim¬ 


ulation figure shows a small northeastern portion of the 
CME impacting Earth in contrast to Figure 6. 

Figures 7b, 8b, and 9b show the modeled magnetic 
held, radial velocity, density, and temperature prohles at 
Venus, Earth, and Mars plotted as color traces for all 48 
ensemble members. Compared with ensemble A (see Fig¬ 
ures 7a-9a), ensemble B results better match the in-situ 
observations at Earth and Venus, particularly the solar 
wind speed at Earth. The modeled magnetic held is lower 
than observed at Venus and Earth because ENLIL is only 
modeling the effect of the pile-up of the plasma and mag¬ 
netic held draping ahead of the CME and does not in¬ 
clude the CME hux-rope. The modeling results show 
that the probability of CME arrival is 38% for Venus, 
35% for Earth, and 100% for Mars. 

At Venus, the ensemble B mean arrival-time is on 9 
January at 05:20 UT with individual arrival-times rang¬ 
ing from 4-10/ — 6 hours around the mean (indicated by 
the grey bar in 7b), giving a prediction error of —0.4 
hours (slight early prediction). At Earth, the ensemble 
mean arrival-time is on 9 January at 20:22 UT with indi¬ 
vidual arrival-times ranging from +10/ — 7 hours around 
the mean, giving a prediction error of +0.7 hours (slight 
late prediction). The standard deviation of the overall 
Kp forecast probability distribution for ensemble B is 
cr = 0.8, and Figure 10b shows that ~ 94% of the predic¬ 
tions fall between Kp = 3 to 5. Either using the average 
Kp prediction of 3 or the most likely Kp prediction of 4, 
the Kp prediction error is 0 or 1. 

At Mars, the ensemble mean arrival-time is on 11 Jan- 
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(o) WSA-ENLIL+Cone Ensemble A Profiles and ACE, WIND Observations ot Earth 



(b) WSA-ENLIL+Cone Ensemble 6 Profiles and ACE, WIND Observatians at Earth 
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Figure 8. Modeled magnetic field, radial velocity, density, and temperature profiles at Earth for (a) ensemble A and (b) ensemble B 
plotted as color traces (color coded by CME input speed) for all 48 ensemble members. In-situ LI observations from ACE are plotted in 
black (red for Bz) and Wind in grey. The vertical dashed red line indicates the CME arrival observed on 9 January at 19:39 UT as a weak 
interplanetary shock. The range of modeled arrival-times for hits is indicated by the grey bar on the time axis. 


Table 2 

Summary of ensemble modeling results for the 7 January 2014 CME at Venus, Earth, and Mars. Column 1: ensemble label. Columns 2-7: median 
CME parameters of speed u, latitude, longitude (HEEQ), major axis half-width a/2, minor axis half-width /9/2, and tilt. The remaining columns 
list the ensemble model mean and spread in arrival-time and the prediction error Aterr = ^predicted ~ ^observed fo^^ Venus, Earth, and Mars. 
Ensemble A measurements and simulations were performed in real-time with a spherical CME shape. Ensemble B measurements and simulations 
were obtained with science level data and use an ellipsoid CME shape with a tilt from CCS measurements. Ensemble B’ is identical to ensemble B 
except a spherical shape is used. Ensemble C half-widths are reduced by 20° compared with ensemble A. Ensemble D longitudes are increased by 

10° compared with ensemble A. 
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uary at 03:13 UT with individual arrival-times ranging 
from -1-26/ — 19 hours around the mean, giving a predic¬ 
tion error of -1-4.7 hours (late prediction). The arrival- 
times at Mars have a large range with most arrivals near 
Ri 22:00 UT on 10 January and some trailing late arrivals. 
For ensemble B, the CME input speed is not the main 
distinguishing factor for arrival-times. Instead, the ear¬ 
lier group of modeled arrivals correspond to CME inputs 
with major half-widths above ~ 40°, and the late group 
have major half-widths in the range of 32°-39°. This 
group of modeled late arrivals suggests that CME major 
half-width is more likely to be above 40°. Ensemble B 


results are summarized in Table 2 and show an overall 
improvement in the CME arrival-time and ifp prediction 
compared with A. 

3.3. Comparing Ensemble Results 

In this section we compare ensemble A modeling results 
with B, which were described above in Sections 3.1-3.2 
and listed in Table 2, and consider other possible reasons 
for ensemble A’s large prediction errors. Often, the CME 
width measurement is poorly constrained, and we note 
that the ellipse major half-widths of ensemble B members 
are on average ss 20° less than the spherical half-widths of 
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(a) WSA-ENLIL+Cone Ensemble A Profiles at Mars 


(b) WSA-ENLIL+Cone Ensemble B Profiles at Mors 
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Figure 9. Modeled magnetic field, radial velocity, density, and temperature profiles at Mars for (a) ensemble A and (b) ensemble B 
plotted as color traces (color coded by CME input speed) for all 48 ensemble members. The vertical dashed red line indicates the observed 
Forbush decrease time. The range of modeled arrival-times for hits is indicated by the grey bar on the time axis. 


ensemble A. In order to better compare ensembles A and 
B, we constructed another ensemble - C ~ such that the 
members have the same CME parameters as ensemble A 
except the half-widths are reduced by 20°. For ensemble 
C the prediction error for the average arrival-time be¬ 
comes —8.5, —11.2, —15.2 hours for Venus, Earth, and 
Mars respectively (compared with —14.9, —19.4, —28.1 
hours prediction error from ensemble A; see Table 2). 
Figure 12a-c show the arrival-time prediction error for 
(a) Venus, (b) Earth, and (c) Mars plotted against the 
CME input speed for all of the ensemble members con¬ 
taining hits. The arrival-time prediction improvement of 
ensemble C (orange triangles) compared with ensemble 
A (black circles) can be seen in this figure. Compared 
with ensemble A, the ensemble C modeled arrival-time at 
Earth becomes later by ~ 4 to 18 hours (8 hours on aver¬ 
age), with the largest differences from CME leg arrivals, 
far from the nose (glancing arrival). The prediction er¬ 
ror improves by « 2 to 15 hours (7 hours on average) 
for Venus and « 8 to 20 hours (13 hours on average) for 
Mars. However, the prediction error improvement is not 
sufficient enough to explain the large early prediction er¬ 
ror of ensemble A solely based on larger widths. 

Next, we consider the possibility of CME deflection or 
the accuracy of the CME longitude measurement. We 
observed that the median CME longitudes of ensemble 
A and B are 20°-30° west of the flare source location 
longitude of 8°. One possible cause for this difference 
was thought to be from the large coronal holes to the 
northeast of the active region (see Figure la) that sub¬ 
sequently produced a high speed stream at Earth on 12 
January reaching Ri850 kms“^. However, Mostl et al. 


(2015) suggest that the CME is channeled due to strong 
nearby active regions to the northeast and open coronal 
field to the west, rather than deflected by the coronal 
holes. The CCS measurements did not show any deflec¬ 
tion in the CME longitude in the measurement range of 
2.1 to 18.5 Rq, which suggests the channeling occurred 
below this height. Similarly, Thalmann et al. (2015) show 
using nonlinear force free modeling that the coronal field 
overlying an active region influences CME speed and di¬ 
rection. Since the CME input parameters making up 
these ensembles are based on coronagraph observations 
subsequent to the lower coronal channeling, and bear¬ 
ing in mind that CMEs in ENLIL are initialized at an 
inner boundary of 21.5 Rq, ensemble A arrival-time pre¬ 
diction errors cannot be explained by this longitudinal 
difference. To illustrate this further, we constructed an¬ 
other ensemble - D - based on ensemble A, such that all of 
the CME input parameters remain unchanged except the 
longitude is increased by 10° for all members. Ensemble 
D is denoted by purple squares in Figure 12 panels a-c. 
The prediction error for the average arrival-time becomes 
—13.3, —17.2, —28.7 hours for Venus, Earth, and Mars 
respectively (see Table 2). Compared with ensemble A, 
the modeled arrival-time at Earth became later by 0.5 to 
6 hours and 0.3 to 4 hours for Venus. For Mars, when 
increasing the longitude, the ensemble D modeled pre¬ 
diction error becomes worse on average by 1 hour com¬ 
pared with ensemble A. Therefore, increasing the longi¬ 
tude leads to a very minor prediction error improvement 
for Earth and Venus and is not nearly sufficient enough 
to explain the large early prediction error of ensemble A. 

The spread in arrival-times from the ensembles shown 
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Figure 10. Probabilistic ifp forecast distribution for (a) ensemble A, (b) ensemble B, (c) ensemble C, (d) ensemble D using ENLIL 
modeled solar-wind quantities at Earth for three IMF clock-angle scenarios 9q = 90° (green), 135° (purple), 180° (orange), and all three 
angles combined 90° - 180° (black) (assuming equal likelihood). Ensemble A shows a 78% chance for Ifp to reach between 6—8 ((t= 1.8), 
ensemble B a 94% chance tor Kp to reach between 3—5 (cr = 0.8), ensemble C an 81% chance for Kp = 3—7 (cr = 1.5) and ensemble D (88% 
chance for Kp = 4-8 (cr = 1.6). 
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Figure 11. Ensemble B global view of the 7 January 2014 CME on 9 January at 18:00 UT: WSA-ENLILJ-Cone radial velocity contour 
plot in the same format as Figure 6 for one ensemble B member predicting a “hit” (CME input parameters: speed of 2061 kms“^, direction 
of 30° longitude, —24° latitude, 49° major half-width, and 31° minor half-width). Panel (d) shows the measured (red) and simulated 
(blue) radial velocity profiles at Earth. Movie available in online version. 
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7 January 2014 CME Ensemble Arrival Time Prediction Error vs. CME Input Speed 



Figure 12. Arrival-time prediction error for (a) Venus, (b) Earth, and (c) Mars plotted against the CME input speed for all ensemble 
members containing hits. Ensemble A measurements and simulations were performed in real-time with a spherical CME shape, and are 
presented by black circles. Ensemble B (green diamonds) measurements and simulations were obtained with science level data and use an 
ellipsoid CME shape with a tilt from CCS measurements. Ensembles C (orange triangles) half-widths are reduced by 20° compared with 
ensemble A. Ensemble D (purple squares) longitudes are increased by 10° compared with ensemble A. 
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in Table 2 depends on a few factors which need to be 
considered when making comparisons. The initial distri¬ 
bution of CME input parameters is the main controlling 
factor. The smaller the spread in the CME parameter 
distribution (smaller error bars), the smaller the spread 
in arrival-times (higher confidence in the arrival-time). 
Ensemble B has larger CME parameter error bars and 
therefore results in larger arrival-time spread compared 
with ensemble A. The spread in arrival-times of glancing 
CME arrivals is expected to be larger (lower confidence 
in the arrival-time) than more central arrivals. In these 
cases, a combination of the CME direction, width, and 
speed can produce a large spread in arrival-times reflect¬ 
ing impacts ranging from near CME center to the legs, 
as seen in ensembles C and D. 

Compared with ensemble A, ensemble B’s CME 
arrival-time and ATp predictions are greatly improved. 
The prediction error from the average ensemble A arrival¬ 
time was —14.9, —19.4, —28.1 hours compared with —0.4, 
-1-0.7, -1-4.7 hours for ensemble B for Venus, Earth, and 
Mars respectively. This can be seen in Figure 12a-c 
where the arrival-time prediction errors at Earth, Venus, 
and Mars are plotted against the CME input speed for 
all members containing hits in ensemble A (black cir¬ 
cles) and B (green diamonds). Comparing Figures 7a-9a 
and Figures 7b-9b the in-situ values and overall mod¬ 
eled arrival-time range is better captured by ensemble 
B. Ensemble A predicted a 78% chance for ifp to reach 
between 6-8, compared with a 94% chance for A"p to 
reach between 3-5 for Ensemble B (see Figures 10a and 
10b). As expected, because ensemble B better captures 
the in-situ observations, the probabilistic Kp prediction 
is closer to the observed value. Ensemble B’s Kp predic¬ 
tions also show improvement compared with ensemble 
D (88% chance for Arp=4-8, ct = 1.6) and ensemble C 
(81% chance for Kp = 3-7, tT = 1.5). We note that the 
most-likely and average predicted Kp from ensemble C 
are similar to ensemble B, however the ensemble C pre¬ 
diction has less confidence (larger standard deviation). 

Next, we consider the background solar wind solution. 
Studies have shown that the CME arrival-times mod¬ 
eled by WSA-ENLIL-bCone are sensitive to the accu¬ 
racy of the background solar wind solution (Lee et al. 
2013; Mays et al. 2015). Temmer et al. (2011) com¬ 
pared WSA-ENLIL modeled background solar wind with 
observations of CME kinematics in interplanetary space 
and found that the ambient solar wind is able to change 
the CME kinematics, for example slow CMEs become 
embedded in the background solar wind. However for 
this period, the model background solar wind predic¬ 
tion prior to CME arrival was close to the observed val¬ 
ues of ~350 kms“^ (see Figure 8). Nevertheless, we 
consider the effect of the GONG daily-updated synop¬ 
tic magnetogram used by WSA to generate the ENLIL 
inner boundary conditions. The original synoptic mag¬ 
netogram used was for Garrington rotation number 2135 
and 40° Carrington longitude on 7 January 2014 at 18:04 
UT. When we change the synoptic magnetogram to a 
Carrington longitude of 50° (7 January at 01:04) and 
30° Carrington longitude (8 January at 14:04 UT) the 
modeled arrival-times for ensemble A do not change by 
more than ±1 hour for Venus and Earth. For Mars, the 
modeled arrivals are shifted earlier by Ri 2 and 6 hours 
for synoptic maps with Carrington longitude 30° and 50° 


respectively. At Mars, the simulations from the different 
maps show a high speed stream on 8 January starting 
earlier in the day and reaching a higher speed, resulting 
in earlier CME arrival-times (« 2 and 6 hours). 

To isolate the effect of the ellipsoid CME shape and tilt 
used in ensemble B compared with the spheroid shape 
used in ensembles A, C, and D, we carried out modeling 
with ensemble B parameters using a spherical cross sec¬ 
tion (and no tilt). This ensemble is labelled B' in Table 
2. Modeled CME arrival-times and Kp predictions for 
this ensemble are very similar to ensemble C (ensemble 
A half-width was decreased by 20°) with a small average 
arrival-time improvement of 1, 3, and 5 hours at Venus, 
Earth, and Mars respectively. Finally, to quantify the 
effect of prior CMEs on the 7 January CME propagation 
we performed additional simulations that included the 
prior 4 January and 6 January CMEs (discussed in Sec¬ 
tion 2). When the prior CMEs are included, the 7 Jan¬ 
uary CME is modeled to arrive ss 1, 2, and 3 hours earlier 
at Venus, Earth, and Mars respectively than when they 
are not included (increasing the prediction error). Taken 
together, these modeling results suggest that the CME 
orientation and directionality with respect to observato¬ 
ries at Venus, Earth, and Mars may play an important 
role in understanding the propagation of this CME. 

4. SUMMARY OF RESULTS 

We presented a series of ensemble modeling results for 
the interplanetary propagation of the 7 January 2014 
CME and arrival at Venus, Earth, and Mars with the 
WSA-ENLIL-bCone model. The 7 January CME was 
chosen for this modeling exercise because its arrival was 
much later than most of the predicted arrival-times on 
the CME Scoreboard, and was far outside the associated 
range of errors for the various methods used. There are 
many potential reasons for the inaccuracy of the WSA- 
ENLIL-bCone prediction; arrival-times depend on the ac¬ 
curate specification of early CME kinematics and geom¬ 
etry, and the model’s ability to reproduce background 
solar wind conditions. For example, a few days later af¬ 
ter the CME arrival (12 January) a high speed stream 
reaching r: 850 kms“^ was observed. Considering the 
large difference (ss 37°) in propagation direction relative 
to the source region longitude, it was thought at first that 
the coronal holes could have deflected the CME from its 
initial path. However, Mostl et al. (2015) suggest that 
the CME is channeled due to strong nearby active regions 
to the northeast and open coronal field to the west rather 
than deflected by the coronal holes. Also, there was a 
delay of more than a day between the in-situ CME and 
high speed stream signatures and the interstitial mea¬ 
surements indicated standard slow wind with no signs of 
compression or interaction. Finally, GCS measurements 
did not show any deflection in the CME longitude from 
2.1 to 18.5 Rq, so the channeling happens in the low 
corona (see e.g. Thalmann et al. (2015)), which indi¬ 
cates that the CME arrival delay is unlikely to be due to 
change in CME longitude in the outer corona. 

The real-time WSA-ENLIL-bCone results from en¬ 
semble A (using default model settings) were system¬ 
atically compared with post-event ensemble simulations 
each testing different effects. We examined the impact 
of CME shape and tilt (ensemble B), CME half-width 
(ensemble C), CME longitude (ensemble D), changing 
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the input magnetogram, and including prior CMEs. The 
prediction error from the average ensemble A arrival¬ 
time was —14.9, —19.4, —28.1 hours compared with —0.4, 
-1-0.7, -1-4.7 hours for ensemble B for Venus, Earth, and 
Mars respectively; ensemble B shows a clear improve¬ 
ment. Because ensemble B better captures the in-situ 
observations, the probabilistic Kp prediction for ensem¬ 
ble B (94% chance for Kp to reach between 3-5) is closer 
to the observed value compared with ensemble A (78% 
chance for Kp to reach between 6-8). ENLIL models 
the effect of the pile-up of the plasma and magnetic field 
draping ahead of the CME and does not include the CME 
flux-rope and a magnetic cloud was not observed in-situ; 
therefore it was not possible to investigate the effect of 
the intrinsic CME magnetic held. However, the different 
effects considered the simulations show that using a tilted 
ellipsoid CME shape with major and minor axis half¬ 
widths would have improved the initial real-time pre¬ 
diction to better rehect the observed in-situ signatures 
at Venus, Earth, and Mars and the geomagnetic storm 
strength. 

The modeling results suggest that CME orientation 
and directionality with respect to observatories at Venus, 
Earth, and Mars play an important role in understand¬ 
ing the propagation of this CME. This suggests that 
CME orientation and shape should be taken into con¬ 
sideration when modeling the propagation of glancing 
CME arrivals, whether due to their asymmetric shape 
in latitude and longitude, or tilt with respect to the so¬ 
lar equator. The 7 January 2014 event illustrates that 
this effect plays a particularly important role for CMEs 
which have a glancing arrival at the location of inter¬ 
est. Many CME propagation models do not currently 
take into account the full 3D geometry of the CME, 
such as tilt, or propagation direction out of the eclip¬ 
tic plane (2D models). For 2D models, the cross section 
of a CCS reconstructed CME along the ecliptic plane 
could be used to improve such modeling results. It is 
important to note that CCS model fits of CMEs are 
not possible without multiple viewpoint imaging such as 
those from STEREO/SEGGHl. Finally, when the mod¬ 
eled CME time of arrival is improved this leads to an im¬ 
provement in the modeled radial velocity of arrival when 
comparing to in-situ observations, which in turn can lead 
to better ENLIL-based Kp predictions. For example, Sa- 
vani et al. (submitted, 2015) show how better estimating 
which part of the CME structure is sampled by an ob¬ 
server during its transit is important for improving Kp 
predictions. 
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